Strong Coupling Electrostatics in the Presence of Dielectric Inhomogeneities 
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We study the strong-coupling (SC) interaction between two like-charged membranes of finite 
thickness embedded in a medium of higher dielectric constant. A generalized SC theory is applied 
along with extensive Monte-Carlo simulations to study the image charge effects induced by multiple 
dielectric discontinuities in this system. These effects lead to strong counterion crowding in the 
central region of the inter-surface space upon increasing the solvent/membrane dielectric mismatch 
and change the membrane interactions from attractive to repulsive at small separations. These 
features agree quantitatively with the SC theory at elevated couplings or dielectric mismatch where 
the correlation hole around counterions is larger than the thickness of the central counterion layer. 

PACS numbers: 87.10.-e, 87.10.Rt, 82.70.-y 



O 

c3 



i 

-a 
c 

o 
o 



> 
in 
in 

00 

cn 

o 

o> 
o 



X 



Biological macromolecules such as DNA, lipid mem- 
branes and proteins are highly charged in water. Elec- 
trostatic interactions play a key role in determining struc- 
ture, phase behavior and specific functioning of these 
macroions in aqueous biological media One particu- 
lar trait of these systems is that their behavior is domi- 
nated, to a large extent, by neutralizing counterions that 
surround them in a diffuse ionic cloud. When present at 
higher valencies, these counterions are known to gener- 
ate strong electrostatic attractions between like-charged 
macroions as observed in numerous experiments and sim- 
ulations 0. These observations stand in stark contrast 
with the traditional mean-field or Poisson-Boltzmann 
(PB) theories, which predict purely repulsive forces [![. 

Consequently, there have been a number of attempts 
to assess corrections to the PB theory using, e.g., cor- 
related density fluctuations around the mean-field distri- 
bution or additional non-electrostatic interactions 0, 0] • 
An alternative approach has been developed recently 
0, 0, 0, 0, 0, HI, which leads to the so-called strong- 
coupling (SC) theory; it is known to become exact in 
the limit of high macroion charge, large counterion va- 
lency, low medium dielectric constant or low tempera- 
ture, where the PB theory breaks down. 

While the strong-coupling phenomena are well- 
understood in the framework of the SC theory, it still 
remains a challenge to predict the behavior of realistic 
biophysical systems, which, among other things, exhibit 
highly inhomogeneous dielectric structure. The large dif- 
ference in static dielectric constants of water (e ~ 80), 
being the most common solvent, and the non-polar moi- 
eties (e ~ 2 — 5), comprising the molecular interiors of 
proteins, lipid membranes and DNA, leads to substantial 
differences in interactions between charges in most com- 
mon biological environments. Moreover, the presence of 
a combination of both aqueous and hydrocarbon regions 



in the immediate surrounding (as, e.g., in the case of two 
or more interacting lipid membranes) leads to a more 
complex pattern of image charges that are induced by 
multiple dielectric discontinuities in the system. 

In this Letter, we consider a system of two like-charged 
membranes of finite thickness in a medium of higher 
dielectric constant and determine the counterion distri- 
bution and the interaction between the membranes by 
means of both Monte-Carlo simulations and a gener- 
alized SC theory. Recent studies show that image in- 
teractions in highly charged systems could result in re- 
markable effects such as in the surface adsorption of 
flexible polyelectrolytes and the charge inversion 
of macroions [HI]. In the slab geometry, the dielectric 
discontinuity effects have been studied in both weak- 
coupling 



9 . [lip. 13l. 14 . 15, [l6[ and strong-coupling 



regimes (9, 3" 16, 17 1. These works however deal with 
one or two semi-infinite slabs and do not consider the 
finite thickness of the membranes or the multitude of im- 
ages produced in the two-slab system. While accounting 
only for the first-order images serves as an accurate ap- 
proximation at low couplings [9L[T3. 14 1 , the same approx- 
imation as we show breaks down on a qualitative level 
at elevated couplings, where the higher-order induced 
images play an essential role. For vanishing membrane 
thickness (no dielectric discontinuity), the long-range SC 
attraction mediated by counterions is dominant and leads 
to a tightly bound state between two like-charged mem- 
branes [1, 0]. This well-known picture changes as the 
membrane thickness and the solvent/membrane dielec- 
tric mismatch increase, leading to enhanced repulsive 
(and even change from attractive to repulsive) interac- 
tion between membranes at small separations. This SC 
repulsion depends strongly on the dielectric mismatch 
contrary to the (repulsive) PB interaction which is not 
affected by the image charges. 
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FIG. 1: Schematic view of two charged membranes and coun- 
terions in between. Here 2a is the closest-approach distance 
(leaving out the counterion diameter from the actual distance) 
and counterion excluded-volume repulsions are neglected. 



Let us consider two charged parallel membranes of 
thickness b and dielectric constant e± with their inner 
surfaces being located at z = ±a and bearing uniform 
charge density —ere (Fig. [I}. The system is immersed in 
a solution of dielectric constant £2 containing +c/-valent 
counterions. Apart from the direct contribution from the 
charged surfaces, the electrostatic potential experienced 
by each counterion also involves contributions from im- 
age charges induced at each solvent/membrane interface. 
We shall consider the contribution from all image charges 
in the present geometry. To this end, we have developed 
a new numerical algorithm that enables one to compute 
efficiently the electrostatic interactions in the presence 
of multiple dielectric discontinuities . Assuming two- 
dimensional lateral periodicity (in the plane of the mem- 
branes), we combine the image charge method with the 
MMM2D summation technique to obtain fast-converging 
series for evaluating Coulombic interactions 17, 

The state of the system may be described in terms 
of the following dimcnsionlcss parameters: the dielectric 
jump parameter A = (£2 — £i)/(£2 + £1), the rescaled 
half-distance a = a/fj,, the rescaled thickness b = b/fj,, 
and the electrostatic coupling parameter 5 = q 2 i~£,lp, 



where fx = l/(2Ttq£Ba) is the so-called Gouy-Chapman 
length and £b = e 2 /(47T£2£ofcBT) is the Bjerrum length 
[H, |2] • The coupling parameter S determines the strength 
of electrostatic correlations: for dielectrically homoge- 
neous systems, the SC behavior dominates typically for 
H > 10, while the PB description is found to be valid at 
small couplings S < 1 0, 0, Q . 

For a typical membrane with o~e ~ l~i i n water (£2 ~ 
80 at room temperature), we have S ~ 3, 25, 85 and 200 
by choosing counterion valency as q = 1, ... ,4, respec- 
tively. Higher couplings may be obtained by using less 
dielectric solvents (e.g., mixtures containing methanol, 
£2 ~ 33) since H ~ Here we shall focus mainly 

on the SC regime with large S and vary A (at fixed 
£2 = 80) in the range A = (no discontinuity, £1 = 80) 
to A = 0.95 (water/hydrocarbon interface, £1 = 2). 

In Fig. [2p,, we show simulated counterion density pro- 
files between two membranes of thickness bj \i = 100 and 
coupling parameter 5 = 100 as A varies (symbols). As 
seen, counterions are strongly depleted from the vicinity 
of the charged surfaces and accumulate around the mid- 



plane as the dielectric mismatch is increased-a trend ob- 
served also in studies that include only the first induced 
images 0, [I2J]. This behavior is intimately connected 
with the electrostatic correlations, i.e., it will be absent 
in the PB limit S — > (see Fig. Intuitively, for 

£2 > £1 counterions have the same sign as their images 
and the counterion-image repulsion provides the mecha- 
nism for the foregoing observation in the SC regime. 

In general, the origin of SC phenomena goes back to 
the fact that at large couplings, S 3> 1, counterions 
strongly repel each other and tend to form a highly cor- 
related quasi-2D layer close to a charged surface. Thus, 
individual ions become increasingly "isolated" as they are 
surrounded by a large correlation hole of size a\ ~ q/a 
(or a±//i ~ vE) as follows from the local electroneu- 
trality condition 0. As a result, the single-particle in- 
teraction between individual counterions (including their 
own images) and the charged membranes, u(z), is ex- 
pected to dominate and lead to a barometric number 
density profile, psc( z ), for counterions; in rescaled units 
p sc (z) ee £sg£ = A(a)e-<*\ where ~z = z/n and 

1/A(a) = J^~dz psc(z). For A = 0, this result indeed 
follows as an exact limiting (S — > 00) result from a sys- 
tematic 1/S-expansion 0, Q, giving a uniform profile 
Psc(z) = ^ (with u(z) — and A(a) = in agree- 
ment with our simulations (horizontal line in Fig. [2^). 
Higher-order terms in general involve corrections due to 
counterion-counterion interactions at finite 5 Q. 

For an inhomogeneous system with A > 0, we gener- 
alize the 1/S-expansion method by taking into account 
the presence of all four dielectric boundaries (Fig. [1]) and 
obtain the leading-order (S — > 00) SC density profile as 



psc(z) 
where Aq = , 
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(1) 



l-cxp(-2Qb) 



This result is plotted in 



l-A 2 exp(-2Qf>) ' 

Fig. [2k (solid lines) and is in excellent agreement with 
our data (symbols) for all dielectric jump parameters at 
S = 100. Qualitatively, the bellshaped profiles may be 
approximated for large b by the first-image expression 



Psc(z) - A(a)e~ 
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The SC pressure acting on each membrane may be cal- 
culated from Psc = Psc / (2tt£bo- 2 kBT) = —dFsc/da, by 
evaluating the SC free energy ^sc (per /cbT and number 
of counterions) using standard SC methods 0, [lj| ; hence 



^SC = a - In 



cosh 2Qz + Aq 
"Jo a( * A -l e 2Q 



-2Qa 



Q e— _ (3) 



The first term is the SC energetic attraction induced by 
single counterions between the two surfaces when there 
are no images (A = 0) [8;, 9\. The second term includes 
the counterion confinement entropy as well as the image 
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FIG. 2: a) Rescaled counterion density profile between two like-charged dielectric membranes for H = 100, b/fi — 100, a/fj. = 1 
and as the dielectric mismatch A is increased (from bottom to top at z — 0) . b) Rescaled inter-membrane pressure as a function 
of A for b/fi — 100, a/fj, = 1 and as S is increased (bottom to top); inset shows the details in the logarithmic scale, c) The 
mid-plane (z = 0) counterion density, p(0), as a function of the coupling parameter 3 for a/fj, = 7, b/fi = 100 and as A is 
increased (bottom to top). In a)-c) lines are the SC predictions (see the text) and symbols are simulation data. 



charge contributions. For A = 0, this second term re- 
duces to the entropic In a term as expected and the SC 
pressure follows as Pqc = —1+1/5, featuring repulsion at 
small separation, a < a*, and attraction at large separa- 
tion, a > a*, with a*//i = 1 being the stable equilibrium 
(zero-pressure) half-distance. Upon increasing A, the 
inter-membrane pressure becomes increasingly more re- 
pulsive at small separations due to the counterion-imagc 
repulsions; it may be increased by up to an order of mag- 
nitude as shown by the data in Fig. 0) (symbols). This 
behavior is captured quantitatively by the SC pressure, 
Psc (solid lines, obtained from Eq. ©) for high enough 
couplings at all values of the dielectric jump. 

So far, we have considered only small membrane sep- 
arations a ~ /i, where good agreement with the SC pre- 
dictions is achieved for S ~ 100 and larger couplings. 
Formally, the SC theory is exact for all separations when 
H — * oo. But for finite 5 as is the case in simulations, 
this limiting single-particle description is expected to per- 
form poorly at large separations, where the sub-leading 
counterion-counterion interactions play a more signifi- 
cant role [f|. For A = 0, the SC validity regime is iden- 
tified as a <C a± or a <C VS 0, i.e., when the surface 
separation is smaller than the counterion spacing, a±, 
defined above. In the presence of images, counterions are 
crowded in a thin layer of thickness (or the full width 
at half maximum) 6 Z < a around the mid-plane (Fig. 
[5^). From Eq. @ the thickness 5 Z may be estimated 
as 6 Z = a[l - (1 + 2aln2/SA)- 1 ] 1 / 2 . The SC theory is 
expected to hold when a± is larger than 8 Z , i.e., when 
counterions form a quasi-2D layer. Hence, we find the 
generalized criterion 8 Z <C a±_ or (in rescaled units) 

5 3 «a»+A» 2 . (4) 

Thus, for the SC theory to be valid at larger a, a larger 
coupling parameter or dielectric mismatch is needed. 
Moreover, this predicts that the SC theory remains valid 



for a larger range of separations when A is nonzero. 

In Fig. [2fc we show the mid-plane density /5(0) as a 
function of 5 for a larger half-distance a/fi = 7. The 
SC theory (solid lines) is no more valid for 5 < 100, 
largely overestimating the simulated density (symbols). 
This reflects the fact that counterion-counterion repul- 
sions absent in the leading-order theory are not negli- 
gible and tend to drive the counterions away from the 
crowded mid-plane toward the surfaces. The same trend 
also transpires in the pressure data for a/fi = 7 (Fig. [3k), 
where good agreement is obtained only for 5 > 500. 

The theoretical and numerical results are generally 
found to be in better agreement for larger A as expected 
from Eq. (gj); see Figs. [5b (inset), [2k and [3^ (inset). 
The change from attractive to repulsive interaction upon 
increasing A is more clearly seen in Fig. [Sk- 
in the presence of dielectric inhomogeneities not only 
the strength of the inter-membrane repulsion becomes 
larger but also its range increases, pushing the attrac- 
tion regime and the equilibrium separation, a*, to larger 
distances. When A = 0, as noted above, a* = a*//i = 1 
in agreement with our simulations in Fig. [3]d. For A > 0, 
a* shows a monotonic increase with the coupling param- 
eter. While the qualitative trend is reproduced by the SC 
theory, the quantitative agreement between simulated a* 
(symbols) and the theoretical predictions (solid lines) is 
obtained only for small to intermediate A. In general, 
the SC theory underestimates the a* values (by < 20% 
in Fig. [3}d). This discrepancy reflects the interplay be- 
tween counterion-image repulsions (included in both the- 
ory and simulations) and counterion-counterion repul- 
sions (included only in the simulations): as A increases, 
the former grows pushing a* to larger values beyond the 
regime of validity of the SC theory (Eq. ([4])), where the 
latter effects need to be accounted for as well. A similar- 
trend is found when membrane thickness is increased. 

Note that while the weak-coupling results are less sen- 
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FIG. 3: a) Same as Fig. [2b but for a/ [i = 7. Inset shows the details for small pressures, b) Simulated equilibrium half-distance, 
a*, as a function of H for = 2 and as A is increased. Lines are SC predictions (from minimization of Eq. @). c) Simulation 
results for a* as a function of b with all images (circles) and the first-order images only (squares) for S = 50 and A = 0.95 (the 
latter shows a very weak increase with b); here lines are guides to the eye. The surfaces attract for a > a* and repel for a < a,. 



sitive to the dielectric jump effects, a remarkable depen- 
dence on A is found in the SC regime (Figs. [2b-c and 
EK)- Since the SC interactions are long-ranged (within 
the regime defined by Eq. ((4])), a larger number of in- 
duced images arc needed to be taken into account in the 
SC regime. The first-order-image approximation [^, [l2j 
fails qualitatively at high couplings and particularly for 
thick membranes as it always predicts a strong attraction 
and a remarkably different (smaller) shown in Fig. 

|3b. 

In conclusion, we have shown that the generalized SC 
theory can explain the effects of dielectric inhomogeneity 
in the two-slab system at elevated couplings (specified by 
an extended SC criterion), where counterions are strongly 
depleted from the interfacial regions leading to enhanced 
repulsive inter- membrane interactions. Thus in contrast 
to weakly coupled systems, the electrostatics of highly 
coupled systems can be affected significantly by the im- 
age charges. As an experimentally measurable effect, the 
SC images may lead to swelling of the multilamellar lipid 
arrays upon changing the solvent dielectric constant. In 
reality, however, one deals with additional factors such 
as smooth dielectric profiles and discrete surface charges 
14 . TR 17 1, which offer interesting problems for a more 
detailed study in the future. 
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